The fermion mass at next-to-leading order in the HTL effective theory 



M.E. Carrington, A. Gynther 

Department of Physics, Brandon University, Brandon, Manitoba, R7A 6A9 Canada 
and Winnipeg Institute for Theoretical Physics, Winnipeg, Manitoba, Canada 

D. Pickering 

Department of Mathematics, Brandon University, Brandon, Manitoba, R7A 6A9 Canada 

The calculation of the real part of a quasi-particle dispersion relation at next-to-leading order in 
the hard thermal loop effective theory is a very difficult problem. Even though the hard thermal loop 
_ effective theory is almost 20 years old, there is only one next-to-leading order calculation of the real 

\ part of a quasi-particle dispersion relation in the literature In this paper, we calculate the fermion 

mass in QED and QCD at next-to-leading order. For QED the result is M = eT/^/^,{l - (1.427 ± 
0.02)e/47r) and for QCD with N[ = 2 and iVc = 3 we obtain M = gT/^/Q{l + (1.867 ± 0.02)sr/47r). 

. PACS numbers: ll.lO.Wx, 11.15.-q 
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CN ! I. INTRODUCTION 

■ It is well known that the behaviour of an elementary particle becomes modified when the particle propagates 
in a medium. The particles become "dressed" by their interaction with the medium, and one speaks of collective 

O ^1 modes, or quasi-particles. One studies these collective modes by looking at the corresponding thermal propagators. 

1^ ' The behaviour of the quasi-particles is deduced from the analytic structure of the propagator. In [3| it was shown 

, from general principles that the singularity structure of certain components of gauge and matter propagators are 
gauge- independent, when all contributions of a given order are systematically taken into account. 

, The calculation of dispersion relations for soft quantities at next-to-leading order (NLO) in the hard thermal loop 

J> ■ (HTL) effective theory is notoriously difficult. There are several calculations of damping rates at NLO. The soft static 

\ gluon damping rate was calculated in Ref. The damping rate of a soft static quark was calculated in Ref. 

■ Calculations of masses and oscillation frequencies require the real part of the dispersion relation, which is considerably 
I more difficult to obtain. There is only one complete calculation in the literature of the real part of a quasi-particle 

■ dispersion relation at NLO: the pure glue plasma frequency in the long wavelength Hmit was calculated by Schulz [ll]. 
VO ' In this paper we calculate the fermion mass, in QED and QCD, at NLO. It is straightforward to obtain the result for 

: QCD from the corresponding result for QED by adjusting the HTL masses and including an overall factor Cf in the 
quark self energy. 

For soft static electrons with momenta = (go ^ eT, q — 0), the mass and damping rate of the quasi-particle are 
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obtained from the solution of the equation 



det(^ --EretiQ)) = 0. (1) 

qo=M-ij 

The fermion self-energy can be decomposed in the usual way: 

llret = T^Slel+T-gSl:^*. (2) 

Since we have taken g = 0, the only non-zero component is S^gi(go) = Tr (7°I]ret('5))/4. From now on, we suppress 
the superscript '(0)' to simplify the notation. In addition, we suppress throughout the subscript 'reV indicating the 
retarded component of the self energy. Using this notation we write the dispersion relation as 

9o-S(9o) - 0. (3) 

qa=M-t'f 

At leading order (LO), the self energy is given by the famihar HTL result, 

ReEHTL = -— , ImSHTL^O. (4) 
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Substituting ^ into ^ we obtain the leading order results for the mass and damping rate: 

m,:=Af(°) = ^,7^°'=0. (5) 

We are interested in obtaining NLO corrections to these results. To obtain these NLO corrections, we expand the 
dispersion relation around the lowest order (LO) solution in Eqn. ([5]), keeping contributions to linear order in NLO 
quantities. The resulting equations are particularly simple because of the fact that the imaginary part of the LO HTL 
self energy is zero. The real and imaginary parts of the dispersion relation give 

M(i)-M(i)ReS]ijTL[^^^°']-ReS]jvLo[M*°^] = 0, (6) 

Using ^ and ^ we obtain 

M« = iReI]A,io[M(0)],7(i)=-ilmI]^io[M(0)]. (7) 

To obtain the NLO mass and damping rate from ((?]) we must calculate the NLO self energy. The original paper 
by Braaten and Pisarski identified three potential contributions. They are: (1) corrections to the LO result for 
the 1-loop diagram obtained by expanding to next order in the ratio of the soft external momentum to the hard loop 
momentum; (2) contributions to the 2-loop diagrams from the region of the phase space that corresponds to both 
loops hard; and (3) contributions to 1-loop diagrams with soft loop momentum, and all propagators and vertices 
replaced with HTL effective ones. The power counting arguments of Braaten and Pisarski refer to the maximum 
possible contribution from each type of term. The actual contribution may be lower order for kinematical reasons, or 
because of some cancellation between different integrals. 

The full NLO contribution is contained in the dressed 1-loop diagrams shown in Fig. [H where the dots on the 
vertices indicate the sum of the bare vertex and the HTL vertex. These diagrams contain all of the contributions 
identified by Braaten and Pisarski, if the loop momentum is integrated over the full range from zero to infinity. The 
integral will also contain subleading contributions that are suppressed by powers of coupling. 




FIG. 1: The diagrams that contribute to the self-energy up to NLO. Wavy lines indicate HTL photons and solid lines are HTL 
fermions. The vertices are defined in Eqn. l[22|l . 



It has been demonstrated that the NLO contribution to the dispersion relation from the diagrams in Fig. [T]is gauge 
invariant. This result is obtained by using the fact that the HTL vertices and propagators satisfy the usual Ward 
identities. One finds that the gauge dependent contribution to the NLO fermion self-energy is proportional to an 
integral times the square of the inverse propagator S'^'^^Q) = Q — Shtl(Q), which vanishes on the mass shell. In 
it was pointed out that a straightforward evaluation of the integral produces mass-shell singularities that cancel 
the contributions from the two inverse propagators, and give a finite gauge dependent contribution to the damping 
rate. This problem was resolved by Rebhan [9] who showed that the integral must be regulated before the mass 
shell is approached. Using this procedure one finds that the position of the pole is gauge independent, and the gauge 
dependence occurs only in the unphysical residue. 

The integral that corresponds to the diagrams in Fig. [T] has the general form 

S((Jo) = dpo / dpJ^{qo,po,P,mf,nbipo),nf{po)) ■ (8) 

J J «0=m/ 
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The factor of in front of the integral is the expHcit factor coming from the two vertices. The integrand is obtained 
by combining HTL propagators and HTL vertex functions, and thermal distribution functions (defined in (fT2l) ). The 
HTL propagators and vertices depend on the 4-momenta and the HTL fermion mass m/. As explained above, we 
substitute the LO result qo = nif in order to extract the NLO contribution. 

We begin by noting that, if the integral in Eqn. ^ is dominated by the part of phase space that corresponds to 
p-soft, we can expand the thermal distribution functions and use: ni,{p) — > T/p and nf{p) — > 0. After expanding the 
distribution functions, one can extract a factor of the temperature, and scale all remaining variables by the LO mass 
TO/. The result has the form 



^NLO 



= e^T • I, 



(9) 



where I is a dimensionless integral that can be calculated numerically. 

For the imaginary part of the self energy, we have explicitly calculated the 2-loop contributions and checked that 
the integral in ^ is dominated by the p-soft region of the phase space. The numerical calculation of the integral 
represented in ^ has been done previously 0, The result is 



7QED 
7QCD 



— . (1.35), 



(10) 



^TCf 
An 



(1.41) for = 3, Nf = 2, Cf = 4/3. 



In this paper we calculate the real part of the self energy by evaluating numerically the integrals that correspond 
to the diagrams in Fig. (TJ without expanding the distribution functions. We extract the numerical coefficients of the 
NLO terms by extrapolating to small values of the coupling constant. The result of this computation is 
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1 + (1.867 ±0.02) 



(11) 



{Nf ^2,Nc = 3). 



II. NOTATION 



In this section we define our notation and give the integrals that determine the real part of the fermion self energy 
at NLO. We use 

= 2g^i/, Smi' = diag (1,-1,-1,-1). 
The thermal distribution functions are defined as 

^ ■^0^^'^f(p)^^pr^'^B{p) = l + 2nb{p),NF{p) = l-2nf{p). (12) 

In this paper we are only interested in thermal effects and consequently we ignore zero temperature pieces of the self- 
energy. We use capital letters to denote 4-momenta: K = {ko,k). We take the external momentum to be Q = (907O) 
and the loop momentum is P = {po,p). We write R — P + Q so that we have r — p. Retarded propagators and self 
energies are obtained from Po ^ Po + and advanced functions from Pq ^ pa — *£• In QED, to leading order the 
fermion and photon thermal masses are 

= -—,ml^—^, (13) 



and in QCD 

m? - '^C,. „4 ^ ^ (iV. + \n,) , & . 4-^. (14) 



where Nc is the number of colours and Nf is the number of fiavours. 
The HTL self energy is 



Shtl(Q) = ^Qo(go,9), Qo(go,9) = ^lnf^^), (15) 
q 2 \qo-" ' 



4 



and we define 

SHTL(d)(Q) = 2zImES=^L(Q), Shtl(s)(Q) = 2ReSS^^L(Q)- (16) 
The HTL fermion propagators are written as 

S{R) = l{S+{R)ijo-f) + ls^{R){lo + f)), (17) 

^+(i?0 = 



2r (m2 + r (r - tq)) + In (7^) (r - tq) ' 



2r 



2r (m) + r{r + tq)) - In (r + m) 

We use the covariant gauge and write the photon propagator in terms of transverse and longtitudinal components 
(recall that p = r), 

D^P) = Pi^MP)+P^,^DUP), (18) 

1 P2 I 

Dt{P) = ^yr^, r, D^{P) 



P^-G(po,r)' ' P^ ~ F{po,r)' 

Furthermore, we define the discontinuities and the principle parts as 

d±(R) = 2ilm5'4(i?), V±{R) ^ReS^^tiR), (19) 
dr/LiP) = 2iImi?^/,^(P), VT/L{P)^ReD^/,\ 
Here, each discontinuity contains a pole contribution and a cut contribution: 

dT/L{K) = -2711 ^ nZT/L{'^T/L{k),k) 6{ko - nujT/L{k)) - 2TTif3T/L{kQ,k), (20) 

ri=±l 

rf±(^) = -27riZ{ko,k) [d{ka-LU±{k))+S{ka+uj^{k))]~2Tri(3±{ko,k). 

Expressions for the functions {/3t, Pl P+, P-}, and the equations from which {lot, lol 1^+, ^-}, are obtained can 
be found using Eqns. (fTZl) and ifTSj) . They are also given in the appendix of [6]. 



III. INTEGRANDS 



At zero temperature, the integral corresponding to the diagrams in Fig. [T] can be written as 

E(0) = -le" j dP{T^{Q,P + Q)S{P + Q)T,{P + Q,Q)D^''{P)+M^„{Q,P,-P,Q)D^,{P)), (21) 

where j dP :— J dpo J d^p, and iS and —iD^^ correspond to the electron and photon lines respectively. We need 
to obtain the corresponding integral at finite temperature. We work in the Keldysh representation of the real time 
formalism. The method we use to sum over Keldysh indices is described in [l3|. The vertices T and M are defined 
in Eqn. l(22l) where P^in indicates the momentum of an incoming fermion, P^out is the momentum of an outgoing 
fermion, and P^j^^ is the momentum of an incoming photon, 

(P^in ; PTpont) — 111 ), (22) 

Mfj_i, (P^in , -P^in : ^7in J P-ipout ) — ^1 {Pip'm -i Pyin ' Pjim P^pout ) • 
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A complete expression for the integrand is derived in Ref. 0|. A reasonably simple form is obtained by rewriting 
the HTL vertices in terms of the self energies, and rearranging the result. There are several tricks that must be used 
to remove the dependence on the HTL vertices. First, one uses the Kubo-Martin-Schwinger (KMS) conditions for 
3- and 4-point functions to obtain an expression that depends only on retarded vertex functions. A complete list of 
the KMS conditions for 3- and 4-point functions is found in For many terms, the HTL Ward identities can be 
used to replace contractions of HTL vertices with the connecting photon momentum by the appropriate difference of 
HTL self energies. There are some terms for which one must use explicit results for the HTL 3-point vertex functions. 
Fortunately, these expressions have a particularly simple form when one of the fermions is not moving. Using these 
techniques, all components of the vertices can be written as simple functions of the HTL self-energy. These self- 
energies also appear in the denominators of the HTL fermion propagators. The general strategy is to rearrange terms 
in the numerators to cancel as many terms as possible with the corresponding factors in the denominators. Significant 
simplifications occur after combining terms and using the mass shell condition = mj. The imaginary part of the 
resulting expression, which determines the damping rate at NLO, has been evaluated numerically in Refs. [1, @|- 

The calculation of the real part of the self energy, which determines the NLO contribution to the mass, is more 
complicated for several reasons. One problem is that second diagram in Fig. [Dproduces pure real tadpole type contri- 
butions that can be dropped in the calculation of the imaginary part. There are additional numerical complications 
that will be discussed in more detail in Sec. HVl We give below the integrals that need to be calculated to obtain the 
NLO contribution to the real part of the self energy We separate terms that contain different combinations of 
delta functions, principle parts, and thermal factors. In addition, we define the operator 



N 



32^ 



dpo / dr. 



(23) 



which will be factored out of all expressions. The integrals are 

12 



— N r^dL{p)NBipo), (24) 

-^Nr^ {-r + qo + raf dL{P)NB {po)V+{R), 
% 

-\n r^{r + qo + r^f dL{P)NB (po) V-{R), 
% 

-^Nr^{-r + qo + r^f d+{R)NF (ro) Vl{P), 
% 

-\n r^{r + qo + r^f d^{R)NF (ro) Vl{P), 
% 

-\n dT{p)NB (po) (2 (-3roi?2 + Amjro + Sr^qo + 6P^qo) - R^J:htl{s){R)) , 
% 

4 N R^Nf iro)VTiP)^iiTL{d){R), 
% 

TV (po - rf (r +P0 + 2qof dT{P)NB (po) V+{R), 

% 

TV (r +po)' {-r + P0 + 2go)' dT{P)NB (po) V-{R), 

% 

TV (po - rf (r + Po + 2<7o)' d+(i?)A^F (^o) Vt{P). 

% 

- ^ TV (r + Po)' (-r + Po + 2q^f d_ {R)Nf (ro) Vt{P). 
% 

The sixth equation in l|24p contains the part of the lowest order result that comes from the Bose-Einstein distribution 
function and the last two equations contain the part of the lowest order result that comes from the Fermi-Dirac 
distribution function. The results for these two terms are 



ReS(di,TVB 
ReI](dL,7'+,TVB 
ReI](dL,7'-,TVB 
ReE(rf+,'PL,TVi. 
ReI](<i_,PL,TVi. 
ReS(dT, A^ij 
ReI](7'T,TVi. 
ReE(dT,n,AfB 
ReI](rfT,P-, A^B 
ReS(d+,7'T, A^i. 
ReS((i_,7'T,TVi. 



Re(dT,AfB) 
Re(d+,7'T, A^f) + Re(d-,PT,TVj.) 



12 
24 



(25) 



where the dots indicate the NLO contribution. 
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IV. NUMERICAL ANALYSIS 



In this section we briefly describe the numeric methods used to get the results given in Eq. (fTTj) . The dimensionless 
integrals to be evaluated numerically, given in Eq. l(24|) , are of the form 



1(9) 



dp / dpo f{g,po,p) , 



(26) 



where g is the coupling constant. The integrand f{g,Po,p) diverges along some curve po — h{p) like [po — h{p)] ^ and 
therefore the integrals must be deflned using a principle value prescription, 



I(q) — lim 



dp 



dpofig,Po,p)+ I dpof{g,po,p) 

lh{p)+e . 



(27) 



The curve h{p) can be computed numerically for each integrand to a high and controllable accuracy. We write 
{h{p) — h(p))/h(p) ^ (5, where h{p) is the numeric estimate of h{p) and (5 is a measure of the error. A numeric 
approximation of the required integrals is obtained as 



dp 



dpo f[g,Po,p) + / dpnj{g,po,p)]+0[e). 

h{p)+e 



(28) 



The integral must be evaluated for a number of different values of e and then extrapolated to e ^ 0. However, for 
numeric stability, one must require S h{p). As a consequence, we have to estimate h{p) to very high accuracy in 
order to get a reliable extrapolation of the limit e ^ 0. This extrapolation is illustrated in Fig. [2] for one value of the 
coupling constant. 
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FIG. 2: Extrapolation of the numeric estimate of the sum of integrals given in Eq. [24]to e ^ as defined in the text (for QCD 
at g = 0.005). 

In order to extract the NLO correction to the thermal fermion mass we must further extrapolate the result to small 
values of the coupling constant. Using the LO result to set the dimensions, we use the ansatz 

^l + a,f+aJ^y (In 1/g + 4) + 0{g% (29) 

where we assume that all coefficients Oj are of order one. Our goal is to determine the coefficient ai. The accuracy 
of the result depends on two things: the accuracy of the numeric estimate of the integrals, and the size of the error 
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that is made by neglecting higher order terms from the ansatz. If we drop terms of order g^, the error in ai can be 
estimated as 



Aai 



27r AS 


+ 


9 ruf 





«2|^(Ml/.9)+«2) 



(30) 



where AS is the error in S(m/) from the numeric estimation of the integrals. Minimizing the error in ai determines 
the optimal range of values of the coupling constant at which the integrals should be computed. We are able to obtain 
an accuracy of the order of AE/m/ « 10~^ which means that using 0.001 < .g ^ 0.006 we obtain an error in the result 
for ai of order Aoi « 0.02 (see Fig.^. 
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FIG. 3: The error in ai as a function of the coupling constant. 



In Figs, m and Owe have plotted the numeric results for I](m/)/m/ and (I](m/) — m/)/m/ • 27r/(7 — ai+0{g) along 
with the best fit (least-square) curve. The best fit results for ai are 



I?™ ^ -1.427 ±0.02, 



(31) 
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1.867 ±0.02 



{Nf = 2,Nc = 3). 
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FIG. 4: The numeric results for E(m/)/m/ along with the best fit curve. 
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FIG. 5: The numeric results for 27r/(?(E(m/) — m/)/m/. The dots are the result of the numeric computations, the solid lines 
correspond to the numeric estimate for ai and the gray regions correspond to the estimated errors as defined in Eq. ll30ll . 



DISCUSSION AND CONCLUSIONS 



The calculation of dispersion relations for soft quantities at next-to-leading order in the hard thermal loop effective 
theory is extremely difficult. Real quantities are particularly hard to obtain. Our results for the next-to-leading 
fermion mass in QED and QCD are given in Eqn (fTTj) . 

It was pointed out in that the subleading correction to the 1-loop HTL contribution (calculated by expanding 
to next order in the ratio of the soft external momentum to the hard loop momentum) is gauge dependent and of 
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order - e3ln(l/e)T. There is a contribution of the same order from 2-loop diagrams where one loop momentum 
is of order T and the other loop momentum contributes a log term coming from an integral of the form / dp^ 
ln(phard/Psoft) ~ ln(l/e). In it was shown that the sum of these two contributions is gauge independent (in the 
class of covariant gauges). Both of these contributions contribute to the coefficient 02 in Eqn. l(29l) and are formally 
included in the integrals given in this paper, which correspond to the diagrams in Fig. [TJ The gauge independence of 
the result in suggests that the full contribution at order 5^ In 1/g might be contained in these diagrams. 
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